Reconsideration of Second Harmonic Generation from neat Air/ Water Interface: 
Broken of Kleinman Symmetry from Dipolar Contribution 
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It has been generally accepted that there are significant quadrupolar and bulk contributions to the 
second harmonic generation (SHG) reflected from the neat air/water interface, as well as common 
liquid interfaces. Because there has been no general methodology to determine the quadrupolar and 
bulk contributions to the SHG signal from a liquid interface, this conclusion was reached based on the 
following two experimental phenomena. Namely, the broken of the macroscopic Kleinman symmetry, 
and the significant temperature dependence of the SHG signal from the neat air/water interface. 
However, because sum frequency generation vibrational spectroscopy (SFG-VS) measurement of the 
neat air/water interface observed no apparent temperature dependence, the temperature dependence 
in the SHG measurement has been reexamined and proven to be an experimental artifact. Here we 
present a complete microscopic analysis of the susceptibility tensors of the air/water interface, and 
show that dipolar contribution alone can be used to address the issue of broken of the macroscopic 
Kleinman symmetry at the neat air/water interface. Using this analysis, the orientation of the water 
molecules at the interface can be obtained, and it is consistent with the measurement from SFG-VS. 
Therefore, the key rationales to conclude significantly quadrupolar and bulk contributions to the 
SHG signal of the neat air/water interface can no longer be considered as valid as before. This 
new understanding of the air/water interface can shed light on our understanding of the nonlinear 
optical responses from other molecular interfaces as well. 
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I. INTRODUCTION 

Since Shen and his co-workers established that the 
second harmonic generation (SHG) from an interface 
can be dominated with the electric dipolar contribu- 
tion of the interfacial molecular layers, SHG has been 
widely used as the spectroscopic probe for molecular ori- 
entation, structure, spectroscopy and dynamics of liq- 
uid/vapor interfaces as well as other interfaces, because 
SHG is electric-dipole forbidden in the centrosymmet- 
ric bulk media, 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 In the early days of 
the development of SHG for interface studies, it has 
been fully realized that, in applying SHG to real prob- 
lems, it is important to know how to assess the rel- 
ative magnitudes of the interface (local) contribution 
and the still possibly significant bulk electric-quadrupole 
and magnetic-dipole (nonlocal) contributions to the SHG 
signali 5 * 1 ^^^ Even though there were exchanges of de- 
bates on whether SHG can be an effective probe for 
isotropic liquid interfaces, i.e. whether the bulk contri- 
bution is negligible (iSiiLi^iiSiSSi it has been generally ac- 
cepted that it is impossible to separate the bulk contribu- 
tion from the total SHG signalii£ii2j2i Over the past two 
decades, theoretical treatment has shown that there is no 
general solution to this problem, and it is often not known 
a priori in interface SHG and sum frequency generation 
vibrational spectroscopy (SFG-VS) studies whether in an 
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interface system interfacial contribution is dominant over 
that of the bulk, or not. 22,23 Recently, Morita proposed 
a formulation towards computation of quadrupolar con- 
tribution in SFG-VS, which can also be used for SHG. 24 
However, no computational result has been reported so 
far. 

On the other hand, a number of experimental stud- 
ies have indicated that '...in the case of insulators and 
liquids with €2 ~ 3, the bulk contribution may still be 
appreciable.'^ and 'For small adsorbed molecules or 
symmetric molecules, or molecules with no delocalized 
electrons, the electric-quadrupole nonlinearity could be 
dominant if both ui and 2oj are away from resonance 
with electric-dipole-allowed transitions.— Thus, it gen- 
erally nullified the possibility for using non-resonant SHG 
to obtain detailed molecular information for pure liquid 
interfaces, and the SHG study on the neat air/water in- 
terface is one of the benchmark examples. This is also 
why the majority of literatures on SHG applications are 
mainly on near resonant chromophores at various inter- 
faces. 

The first neat liquid interface studied with SHG experi- 
ment is the neat air/ water interface by Goh et aim In this 
and the following works, two basic phenomena of the neat 
air/water interface were observed. One is the significant 
temperature dependence of the SHG signal (the SHG in- 
tensity dropped more than 50% when temperature was 
changed from 283K to 353K), and the other is the bro- 
ken of the Kleinman symmetry, i.e. Xxzx ~ 2-3Xzxx at 
293K.Mi Since it was generally believed that the Klein- 
man symmetry rule for dipolar susceptibility states that 
in the static limit, i.e. far from electronic resonance, 
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metry for the neat air/water interface was believed to 
indicate a significant quadrupole (bulk) contribution. 3 4 
In addition to the broken of Kleinman symmetry, the 
strong temperature dependence was also believed to in- 
dicate significant quadrupolar contribution. The temper- 
ature dependence was further used to separate the dipo- 
lar (interface or local) contributions and the quadrupo- 
lar (bulk or nonlocal) contributions, because the dipolar 
part should strongly depend on the absolute orientation 
of the water molecules, whereas the quadrupolar con- 
tribution should only be weakly dependent on suchA* 
Using a simple model, Goh et al. further estimated that 
the energy that causes a net absolute orientation of wa- 
ter molecules at the interface to be about kkT at room 
temperature^ 

However, the temperature dependence result in SHG 
measurement was questioned by Du et al. with the first 
SFG-VS investigation of the neat air/water interface, 
where the vibrational spectra of the interfacial free OH 
group was directly measured^ Temperature dependence 
of the SFG spectrum of the air/water interface showed 
no significant change over the same temperature range 
from 283K to 353K— As reported, Du et al. also re- 
peated the SHG temperature dependent experiment of 
the neat air/water interface from 283K to 318K, and they 
observed no change of Xxzx within the experimental un- 
certainty of ±10%, while a 50% decrease was expected 
from the Goh et aVs results i&ft^Si 

On the request by Professor Y. R. Shen, Professor K. 
B. Eisenthal's group repeated the SHG temperature ex- 
periment of the neat air/ water interface in January of 
1993, and found that in the same temperature range, 
Xxzx decreases only by less than 10% and Xzxx decreases 
by about 10%, as reported in Du et al.M One of us 
(Hong-fei Wang) was then the student in the Eisen- 
thal group who was responsible for repeating the SHG 
temperature dependence measurement. It was found 
that the strong SHG temperature dependence reported 
previously 3 ** was an experimental artifact due to the un- 
controlled condensation of water vapor on the internal 
surfaces of the two quartz windows of the enclosed tem- 
perature controlled cell made of brass. The condensa- 
tion became visible when the internal temperature of the 
cell is higher than the room temperature controlled at 
22.0±1.0°C, and more condensation was formed when the 
temperature difference became larger before it reached a 
steady state. Thus, the SHG signal drop became more 
pronounced with the increase of the degree of condensa- 
tion until it levelled off around 310K. Heating of the ex- 
ternal surface of the temperature cell with slightly higher 
external temperature than the internal temperature can 
prevent condensation on the internal surface of the quartz 
windows. After getting rid of the condensation, no SHG 
signal temperature dependence was observed beyond the 
experimental uncertainties, which were 7% for the Xxzx 
and 10% for the Xzxx, which is significantly smaller £L 
Thus, the strong SHG temperature dependence reported 
previously** was an experimental artifact, and was not 



consistent with the SFG-VS measurement — Recently, 
with the knowledge that Goh et al.'s SHG temperature 
dependence results were not repeatable, Fordyce et al. 
reported a new measurement that there was a 20% of 
the SHG signal temperature drop of the neat air/water 
interface from 293K to 343K, 12 Careful examination of 
this report makes us believe that without effective con- 
densation control for the temperature cell windows, this 
work might have also been subjected to the same exper- 
imental artifact as that of Goh et al.. Nevertheless, all 
the SHG experimental results confirmed the validity of 
the broken of Kleinman symmetry for the neat air/water 
interface,**.** 2 . 

Because the magnitude of the temperature dependence 
of the SHG signal was believed to indicate the mag- 
nitude of the dipolar (local) contribution to the SHG 
signal, 3,4,12 no apparent temperature dependence could 
lead to the conclusion that the apparent dipolar con- 
tribution to the SHG signal from the neat air/water 
interface is negligible. In addition, the broken of the 
Kleinman symmetry was also an indication of significant 
quadrupolar (nonlocal) contribution to the SHG signal 
from the neat air/water interface. Therefore, it seems 
reasonable to conclude that all the SHG signal from 
the neat air/ water interface was quadrupolar contribu- 
tion in nature! Indeed, Bloembergen et al. did conclude 
as early as 1968 that the 'SHG reflection from media 
with inversion symmetry is described rather well by the 
quadrupole-type nonlinear properties calculated for the 
homogeneous bulk material with a abrupt discontinuity 
at the boundary— However, Bloembergen et al.'s dis- 
cussion was on experimental results and theoretical anal- 
ysis of metal and alkali halides surfaces, which possess 
relatively high density of valence electrons, but not on 
surfaces of normal liquid. Such a strong statement on 
metal and alkali halides surfaces by Bloembergen et al. 
was firstly questioned by Brown et al's SHG experiment 
on adsorbed surface layer on silver in 1969;** but only 
fundamentally modified more than a decade later, af- 
ter Y. R. Shcn's group demonstrated systematically in 
the early 1980's that surface SHG can be treated as ra- 
diation from a nonlinear polarization sheet induced by 
an incoming wave at the surface, besides comparable 
quadrupole (bulk) contributionsi 5 i 30 i 31 ' 32 i 33 This devel- 
opment undoubtedly assured the surface monolayer sen- 
sitivity of surface SHG for all kinds of surfacesi 5 i 30 i 31 i 32 i 33 
These works by Y. R. Shen and his colleagues heralded 
the beginning of systematic interface studies with nonlin- 
ear optical techniques)" 5 including studies of the longtime 
nail-biting liquid interfaces** These twists in the histori- 
cal development of surface SHG demonstrated the com- 
plexity and difficulty on the evaluation of relative bulk 
and surface contributions to surface SHG signal. 

Now, if the above conjecture on the pure quadrupolar 
(bulk) contribution to SHG from neat air/water inter- 
face could be held true, the observed SHG signal could 
have no surface specificity and should not have been used 
to study the molecular orientation and structure of the 
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air/water interface. And this can also be the same for 
other liquid interfaces. The consequence is that the whole 
effort of using SHG for non-resonant liquid interface stud- 
ies should have been nullified. Indeed, in the past decade, 
using SHG to study neat air/liquid and liquid/liquid in- 
terfaces has been ra.re£^H This situation would persist 
unless the widely held assumptionjiiLiS that the broken 
of the Kleinman symmetry is an indication of significant 
quadrupolar (nonlocal) contribution to the SHG signal 
from an interface, is generally not true. Examination of 
the history of the Kleinman symmetry of the nonlinear 
microscopic polarizability and macroscopic susceptibility 
tensors indeed raised more questions on the general va- 
lidity of this symmetry relationship which first discussed 
by Kleinman in 1962i2 5 - 3 ii 

The validity of the Kleinman symmetry of the sec- 
ond order nonlinear susceptibility tensors was thoroughly 
discussed by Franken and Ward as early as in 1963, 35 
and most recently by Dailey et aZ. 36 As pointed out by 
Franken and Ward, Kleinman suggested that for a mate- 
rial transparent (far from resonance) to the fundamental 
and the SH frequencies, an additional symmetry condi- 
tion of Xijk = Xjik exists, where the first index is for 
the SH frequency, and the second and third indexes are 
for the fundamental frequency^ 5 - To translate this for 
water, Xxzx = Xzxx should hold be true if the Klein- 
man symmetry is valid when both fundamental and the 
SH frequencies are far from resonance with the water 
molecule transition frequencies. It was well recognized 
by Franken and Ward that the presence of dispersion 
can destroy the validity of the Kleinman symmetry, and 
using a quantum mechanical formulation of the dipolar 
susceptibility tensors they further pointed out that even a 
few percent dispersion can seriously damage the validity 
of Kleinman symmetry. It was further demonstrated by 
Franken and Ward that for the case of quartz with the 
fundamental frequency at 694. 3nm of the Ruby laser, 
even though the dispersion is only about 2%, the bro- 
ken of the Kleinman symmetry can be as significant as 
30% by considering the dipolar contribution onlyi^ 5 . This 
undoubtedly suggested a general breakdown of the Klein- 
man symmetry for general materials and molecules. This 
is also why Franken and Ward^ as well as Ron ShenjSI 
called it as Kleinman's conjecture instead of Kleinman 
Symmetry^ Therefore, in the study of real molecular 
systems, Kleinman symmetry should generally fail un- 
der the dipole approximationiSiiS Recently, using the 
Kramers-Kronig dispersion relations, Dailey et al. again 
discussed this general failure of Kleinman symmetry rela- 
tionship in nonlinear susceptibility tensors of molecular 
systems, and they further listed numerous experimen- 
tal disagreements with the Kleinman symmetry relation- 
ship under non- resonant conditions in the literatures^ 
Thus, it is clear that the broken of Kleinman symme- 
try can never be used as the criterion for the existence 
of quadrupolar (bulk) contribution to SHG signal. How- 
ever, the validity of the Kleinman symmetry had been 
overvalued^^iiS due to the lack of awareness about the 



original work by Franken and Ward. 35 

In Franken and Ward's discussion of quartz crystal, 
the macroscopic symmetry of the crystal is the same 
as the microscopic symmetry of the unit cell. How- 
ever, for molecular system, the macroscopic symmetry 
of the molecular system and the microscopic symme- 
try of the molecule are usually not identical, especially 
for the liquid interfaces which usually have rotationally 
isotropic symmetry around the interface normal, no mat- 
ter what symmetry each individual molecule may pos- 
sess. Therefore, the connection between the microscopic 
and macroscopic Kleinman symmetry relationship is not 
as straightforward as that for crystals. Recently, Simp- 
son et al. have tried to bridge this gap by working on 
the molecular and surface hyperpolarizability of oriented 
chromophores of different symmetries and different elec- 
tronic transitions^Si^2ii£ These works can certainly pro- 
vide good references for dealing with complicated rela- 
tionships of susceptibility tensors in SHG and SFG stud- 
ies of molecular interfaces. 

For the neat air/water interface, Goh et al. pointed 
out that derivation of the expressions of biaxial molecule 
{C% v symmetry) might help interpretation of the SHG 
signal from the air/water interfaced This was indeed pur- 
sued briefly for the air/water interface^ but it was ul- 
timately carried away by the observed temperature de- 
pendence phenomenon, which seemingly indicated sig- 
nificant quadrupolar contributions to the SHG signal. 3,4 
Besides, the major obstacle in the advance of SHG tech- 
nique for surface studies has always been the lack of mi- 
croscopic understanding of the surface nonlinear suscep- 
tibility Xijfc- 

In this report, after brief discussion of the basic the- 
ory of SHG, we shall first present the discussion on the 
connection between the microscopic polarizability ten- 
sors of the water molecule and the macroscopic suscepti- 
bility tensors. Then we shall discuss the general con- 
siderations regarding the interface and bulk contribu- 
tions to the SHG from an interface. Finally we shall 
present the treatment on the SHG measurement with 
the proper formulation of the microscopic polarizability 
tensors of the water molecule, which gives well descrip- 
tion with the dipolar nature of the SHG response from 
the neat air/water interface. The major conclusion is 
that with explicit treatment of the molecular symmetry, 
dipolar (interface) contribution alone can provide satis- 
factory explanation to the observed SHG signal from the 
neat air/water interface. We shall show that not only 
detailed information on molecular orientation can be ob- 
tained from polarization SHG measurement of the neat 
air/water interface, but also such treatment can be gen- 
erally applied to other molecular interfaces as well. This 
new understanding of the air/ water interface can shed 
light on our understandings of the nonlinear optical re- 
sponses from other molecular interfaces. 
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II. THEORETICAL BACKGROUND 
A. Basic Theory 

The basic theory of SHG as a general surface analytical 
probe has been well described in the literaturesi 5 * 4 ^ 3 ^^ 
However, detailed treatment of the experimental data in 
different polarizations has not been usual4244i2£4£i2L2& 

Generally, the SHG Intensity I(2u>) reflected from an 
interface is given below 



I(2w) 

Xeff 
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(1) 



[L(2w) : e(2aj)]. X ■ [L(w) : e(w)].[L(w) : e(w)] 

(2) 



In Equation^ I w is the incoming laser intensity, cq is the 
speed of the light in the vacuum, and £1 is the incident 
angle from the surface normal. In Equation HI x is the 
macroscopic second-order susceptibility tensor, which has 
3x3x3 = 27 elements; e(2cu) and e{u>) are the unit vector 
of the electric field at 2u> and w; L(2cj) and L(w) are the 
tensorial Fresnel factors for 2uj and ui, respectively. 

It is important to realize that Xef f contains all molec- 
ular information of SHG measurement. Actually Xeff 
is the point product of the observation vector \L,(2u>) : 
e(2uj)} and the second order polarization vector — 
X : [L(w) : e(w)].[L(w) : e(uj)]. In the SHG experi- 
ment, the field vectors of the incoming and out-going 
light beams are controlled by the experimenter; once 
these field vectors are fixed, the tensorial Fresnel factors 
are also fixed quantities. So the SHG experiment mea- 
sures the SHG intensity and provides the information of 
the macroscopic susceptibility of the molecular system. 

In SHG experiment, there are three independent polar- 
ization measurement, namely, s-in/p-out (Isp), 45°-in/s- 
out (/450s) and the p-in/p-out (I pp ). Here, in the ex- 
perimental coordinate system (x,y,z), z is the interface 
normal, and we choose xz plane as the incident plane. 
Subsequently, p polarization is defined as polarization 
within the xz plane, and s is perpendicular to the xz 
plane. If the microscopic local field factors of the molec- 
ular layer are considered, the Xeff of these three polar- 
izations should take the following forms 4242^ 



Xef f,sp 
X e fjAb°s 



Xeff 



I'll 



L zz {2u)L 2 yy {uj)l zz {2uj)l 2 yy {u})sv[Y ttXzyy 
L yv {2u))L zz {u)L yy (u) 

Xlyy(2U))l zz (uj)lyy(uj) SUlflXyzy 

L zz {2lo)L 2 xx {uj)I zz {2uj)II x {u) sin ft cos 2 Qxz 
-2L xx (2u)L zz (u))L xx (w) 
xl xx (2u>)l zz (u>)l xx (uj) sinficos 2 ttxxzx 
+L zz (2uj)L 2 zz (lu)I zz (2uj)1 2 zz (uj) sin 3 n Xzzz 



where l xx (tOi), l yy (uji) and l zz (u>i) are the microscopic 
local field factors at frequency ^-, 43 i 50 i 51 The detailed 
expression for the La and la factors can be found 
elsewhere, 43 i 49 i 50 i 51 

Now the intensity ratio Ieff^s/Ieff.sp should be the 
square of the ratio XeffA5°s/Xeff,sp- 



XeffA5°s 
Xeff,sp 



L V y(2uj)L zz (u)) l yy (2u)l zz (u>) x^ 

L ZZ (2Uj)Lyy(uj) l ZZ {2uj)lyy{uj) % 



(4) 



So the difference of the I e ffA5°s an d I e ff,sp can be either 
from the difference of the Fresnel factors and local field 
factors at u> and 2lu, or from the difference between the 
second order susceptibility tensor Xyzy and Xzyy 

As we know, when both u> and 2oj are far from reso- 
nance with the molecular electronic transitions, the dif- 
ference of the Fresnel factors and local field factors at lo 
and 2ui are essentially negligible ^SiiiSi Therefore, gen- 



erally we have 



t {2uj)L zz (u>) 



1 and 



lyy(2u)l zz (uj) 



1 for 



L zz (2U))Lyy(uj) l zz (2u)lyy(u) 

simple liquid interfaces far from resonance. 12 Usually, the 
linear dispersion of non-resonant materials in the optical 
frequency region is only in the order of 1% to 3%i^ 

The non-zero elements of the macroscopic susceptibil- 
ity tensors are determined by the symmetry property 
of the interface. For an achiral rotationally isotropic 
molecular layer, there are only 7 non-zero elements, i.e. 

Xzzz, Xzxx = Xzyy j Xyzy = Xyyz = Xxzx = Xxxz- If 

Klcinman symmetry holds additionally, we have Xzzz, 

Xzxx — Xzyy — Xyzy — Xyyz — Xxzx = Xxxz^^^ m There- 
fore, if the ratio Xeff,45°s/Xeff,sp is significantly different 
from unity for a simple liquid interface, i.e. Xyzy 7^ Xzyy 
from Eq^l the macroscopic Kleinman symmetry must be 
broken. 



B. Microscopic and Macroscopic Kleinman 
Symmetry 

Here we discuss the direct connection between the 
Kleinman symmetry of the microscopic (molecular) po- 
larizability tensor elements (3i mn and the macroscopic 
susceptibility tensor elements Xijk of the molecular sys- 
tem. Here we only consider the contribution from the 
dipolar terms. 

Under the dipolar approximation, it is known that Xijk 
is generally related to (3i mn through Euler transformation 
operation, which is the product of three rotational oper- 
ations of the molecular coordinates systemi^Li&iji 



Xijk 



N 



^ {RilRjmRkn) Pl 



(5) 



//'(• 



(3) 



where N is the interface molecule density, the operator () 
denotes the orientational ensemble average over the Eu- 
ler rotation matrix transformation element R\y from the 
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molecular coordinate A' (a, b, c) to the laboratory coordi- 
nate A (a;, y, z)i 44 i 55 The subscript k) of the Xijk cor- 
responds to the laboratory coordinate (x,y,z), and the 
subscript (I, m, n) of the /3; rml corresponds to the molecu- 
lar coordinate (a, b, c). This expression can also be writ- 
ten into the summation over [3i m n of each molecule in the 
molecular system which contributes to Xijk- 



N 

Xijk — J]] 

1=1 It. 



'Iran 



(6) 



a be 



In the most general case there are 27 second order sus- 
ceptibility tensor elements for Xijk as well as (3i mn - Be- 
cause in SHG the two fundamental frequencies are equiv- 
alent, the relationship Xijk = Xikj and (3i mn = fiinm 
are always valid. So in the most general case there 
are only 18 independent elements for Xijk and f3i mn , re- 
spectively Kleinman symmetry states that Xijk — Xjik 
for the macroscopic susceptibility tensor elements, and 
Pimn = Pmin for the microscopic polarizability tensor el- 
ements. Kleinman symmetry can further reduce the in- 
dependent elements to the number of 10. Of course, the 
number of non-zero tensor elements is much smaller than 
this number. 

Because the microscopic Kleinman symmetry relation- 
ship requires /3i mn = fimim and because the product 
of the matrix elements always satisfy RfiRj m Rf. n = 
RjmRii-R-ln' we immediately have 



JV 



Xijk 



Hran 



2_j RilRjmRknflh 
1=1 lmn=abc 
N 

^ ] ^ t RjmRilRknPrnln — Xjik 



1=1 I 



mn=abc 



(7) 

This simple transformation proves that the Kleinman 
symmetry is invariant under arbitrary rotational trans- 
formation. Because the rotational transformation is lin- 
early reversible, if the macroscopic Kleinman symmetry 
holds, the microscopic Kleinman symmetry should also 
hold, as long as all the molecules in the system can be 
treated identical. This relationship conceptually simpli- 
fies many things when we are trying to compare experi- 
mental results with molecular properties in SHG, as well 
as SFG studies. 

Microscopic Kleinman symmetry is valid for uniax- 
ial rod-like molecule, because such molecule has only 
one non-zero molecular polarizability tensor element (3 CCC . 
Therefore, any macroscopic system consists of uniaxial 
molecule shall observe Kleinman symmetry. It was just 
because the 'ubiquitous' usage of uniaxial molecule model 
in description of nonlinear chromophores, Kleinman sym- 
metry has been generally held true in the past forty 
years, pointed out precisely by Dailey et ai.jSS despite 
the seminal discussion of Franken and Ward in 1963 on 



the general failure of Kleinman symmetry even for quartz 
crystal^ 



C. SHG from Liquid Interface of Water: the 
Molecule with C2v Symmetry 

Here we present the detailed treatment of the SHG 
from neat air/ water interface considering only the dipolar 
contribution. 

The water molecular belongs to Ci v symmetry. There- 
fore, there are seven non-zero microscopic polarizability 
tensor elements in SHG, namely, /3 CCC , (3 caa , (3 aca = /3 aac , 
Pcbb, Pbcb = fibbe fi^^ Here the molecular coordinates 
are defined to have c axis along the Ci axis, and to have 
ac plane as the molecular plane. For the rotationally 
isotropic neat air/water interface, there are 7 non-zero 

elements, i.e. Xzzzi Xzxx Xzyyj Xyzy Xyyz Xxzx 

Xxxz- Because the twisted angle ip and the azimuthal an- 
gle <j) can be considered isotropic, using Eq[S]and average 
over Euler angle ip and (f>, we have 



Xzzz = -y [{flcaa + Pcbb + 2f3 aca + 2[3 bcb ) (COS 9) 

+ {2f3 ccc - f3 caa - /3 cbb - 2/3 aca - 2/? bcb )(cos 3 6)] 
Xzxx = -£-[{20 ccc + (i caa + (3 cbb - 2[3 aca - 2/3 bcb )(cos 9) 

~{2f3 ccc - /3 caa - f3 cbb - 2f3 aca - 2/? bcb )(cos 3 9)} 
N s 



_ S 

Xxzx 7~ 



[(2/3 ccc - f3 caa ~ /3 C 6&)(COS0) 

-(2f3 ccc - (3 caa - /3 cbb - 2(3 aca - 2/? bcb )(cos 3 9)] 

(8) 

Because /3 caa and (3 cbb always go together in these expres- 
sions, if I3 CCC ^ 0, we define r — hb ', and similarly 

we define s = P ac °-+P bcb . Then we have 

Xzzz = N s /3 ccc [{r + 2s){cos9) + (l-r-2s)(cos 3 9)} 
Xzxx = ^(3 ccc [{l + r-2s)(cos9)-(l-r-2 S ){cos 3 9))} 



Xxzx = -y/W(l - r ) ( cos e ) - C 1 - r - 2s ) ( cos3 



And the Xzxx/ Xxzx and Xzzz/Xxzx ratios are 



(9) 



Xz 



XxZX 

"Xzzz 
Xxzx 



{1+r -2s)D - (l-r-2s) 

(1 - r)D - (1 - r - 2s) 
2(r + 2s)D + 2(l-r-2s) 
(1 - r)D - (1 - r - 2s) 



(10) 
(11) 



in which D = (cos 9)/ (cos 3 9) is the orientational pa- 
rameter of the water molecule at the neat air/water 
interfaced It is clear from Eq^J that if r = s, then 
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Xz 



Xxz 



and vice versa. From the definition, 



is equivalent to /3 caa + (3 cbb = (3 aca + (3 bcb , which is the 
necessary condition for Kleinman symmetry. 

When r ^ s, it is clear from Eq^| that the ratio 
Xzxx/Xxzx 1 unless D = 0. However, D = is phys- 
ically impossible, therefore, a trivial solution. 44 Thus, 
Xzxx/Xxzx value not only depends on r and s, but also on 
the value of the orientational parameter D. Only when 
D = 1, we have Xzxx/Xxzx = r/s, which means that the 
macroscopic and microscopic tensor ratios are equal. Be- 
cause the two ratios of Xzxx/Xxzx and Xzzz/Xxzx can be 
obtained from experiment by using Eq|3 and if we know 
the value for r/s, D value can be determined by solving 
the EqUHIand EalTTI 

Here with the special case that when p ccc = 0, only 
the ratio R = r/s = g 222 ^^ is physically meaningful, 
and Eg 1 101 and Eq^] become into the followings. 



Xzxx = (R-2)D + (R+2) 
-R*D + (R + 2) 
2{R + 2)D -2(i? + l 
XxTx ~ -R*D+(R + 2) 



X.XZX 

X z z z 



(12) 
(13) 



It has been known that from quantum mechanical 
treatment, C 2v molecule must have a low- lying excited 
electronic state with a transition dipole moment perpen- 
dicular or close to perpendicular to the molecular dipole, 
i.e. along a or b axis^ This low-lying state belongs to 
the -B-type irreducible representation, and the dominant 
molecular polarizability tensors for the £?-type transition 
can only be either (3 aac = (3 aca (Bi) or (3 bbc = f3 bcb (B 2 ) 
in resonance^ For water molecule, the first two elec- 
tronic excited states belongs to B\ and B 2 symmetry, 
and the transitions are around 140nm-190nm (diffusive) 
and 124nm (strong), respectively^ According to Dai- 
ley et al.fiSi the Kramers-Kronig dispersion relations dic- 
tate that even far from resonance, the largest tensor el- 
ements for water molecule have to be (3 aac = j3 aca or 
flbbc = /3fc c 6, and the Kleinman symmetry has to be bro- 
ken. Importantly, because the first two transitions are 
along a or & directions, /3 CCC is expected to be very small 
or negligible, comparing to f3 aac = (3 aC a, Pbbc = Pbcb, 
(3 caa and p cbb terms, according to the quantum me- 
chanical sum-over-states (SOS) treatment of the molecu- 
lar hyperpolarizabilitiesi^^Si^ Therefore, the broken of 
macroscopic Kleinman symmetry for the neat air/water 
interface reported previously is unavoidable, and can be 
fully formulated using only dipolar contribution as de- 
scribed in this section. 

If Pccc term can be very small or negligible as discussed 
above, Eg 1 121 and Eq^|can be used to extract unique R 
and D values. 



D. Beyond the Dipolar Contribution 

Here we discuss the implication of the above treatment 
on the dipolar contribution to the total SHG signal. 



Detailed treatment on contributions to the total SHG 
signal from the dipole, quadrupole and field gradient 
across the interface layer has been well established in the 
literaturesAiiiiiSiSi The dipolar contribution is usually 
referred as the interface or local contribution, and the 
other two contributions are generally referred as the bulk 
or nonlocal contributions. 14 When the upper phase is air, 
we have, 



total 

Xzzz 



total 

Xzxx 



total 

Xxzx 



Xs,zzz \i Xq ,zzzz) 



e'(2to)e' 2 (uj) 



e{2uj)€ 2 {uo) 
-( + /(^) 

X.s,zxx w X-Q.zzxx) £^2uj^ 

e'(a,) 



(14) 



where xTjk 1 1S the total susceptibility tensor element, 
which includes both the interface and bulk contributions. 
Xs.ijk is the dipolar susceptibility tensor, denoted as Xijk 
in previous sections. 7 represents the contribution from 
the field gradient across the interface. e'(w), e'(2u>) are 
the effective optical dielectric constants of the interface 
layer, and e(uj), e(2tu), are the optical dielectric constants 
of the bulk liquid. Xq ijki 1S ^ ne quadrupolar susceptibil- 
ity tensor components of the bulk medium. 

As we have known, the ratios between the three y'°£ a/ 
in Ea 1141 can be determined directly from the SHG mea- 
surements. However, in principle the relative contribu- 
tions to xljk al from the local and nonlocal terms can not 
be explicitly estimatedii^iSiSi From Eq^]Here we argue 
that if the ratios of the Xijk al terms can be quantitatively 
explained with the dipolar terms, i.e. Xs,ijk terms, the 
contributions from the bulk terms can be very small or 
negligible. At least, based on the discussions we have so 
far, all the previous arguments for them to be comparable 
to the dipolar terms can not be substantiated. 



III. EXPERIMENTAL 

A broadband tunable mode-locked femtosecond 
Ti:Sapphirc laser system (Tsunami, Spectra-Physics) is 
used for reflected-geometry SHG measurement. Its high- 
repetition rate (82MHz) and short pulse-width (80fs) 
make it suitable for detection of weak second harmonic 
signals. Its long term power and pulse- width stabil- 
ity also make it easy for quantitative analysis of the 
SHG dataiS The optical set-up is typical for SHG 
experiment^ The 800nm fundamental laser beam is fo- 
cused on the air/water interface at the incident angle of 
il = 70° , and the SHG signal at 400nm is detected with a 
high gain photo-multiplier tube (R585, Hamamatsu) and 
a photon counter (SR400, Stanford). Typically the dark 
noise level is less than 3 counts /second. The laser power 
is typically 500mw. The efficiency of the detection system 
for p polarization is 1.31 times of that for s polarization. 
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A teflon beaker (diameter around 10cm) is used to con- 
tain the pure water liquid (Millipore 18.2MQ-cm). Op- 
tical polarization controlling, SHG data acquisition are 
programmed and controlled with a PC. It takes about 
300 seconds to collect the polarization curve. The room 
temperature is controlled at 22.0 ± 1.0°C, and the hu- 
midity of the room is controlled around 40%. 



IV. RESULTS AND DISCUSSIONS 
A. Calculation of SHG Intensities 



TABLE I: The Xzxx/Xxzx and Xzzz/Xxzx ratios from different 
research groups. The calculated R, D and 9 values according 
to EqEUand EatTill are also listed. 





Eisentha! 3 


Frej^ 


Girault 9 


This work 


Xzxx 


0.46 


0.48±0.02(?) 


0.41 


0.30±0.03 


Xzzz 


1.20 


3.6T0.5 


2.56 


2.56T0.06 


R 


0.66 


0.81T0.03 


0.74 


0.69T0.02 


D 


1.39 


1.84T0.07 


1.69 


1.72T0.02 


e 


32.0° 


42.5° ±1.2° 


39.8° 


40.3°±0.4° 



Here we present the analysis of the polarization SHG 
data with the dipolar formulations in Section II. 

The SHG signal in both p and s polarizations are mea- 
sured with automatic scanning of the input polarization 
in the full 360°, as presented in Fig^ Because the SHG 
signal from a neat air/liquid interface is usually very 
small, full range scanning of the input polarization can 
help improving the accuracy of the SHG intensities mea- 
surement in the desired sp, 45°s and pp polarizations. 
The SHG signal is especially small for sp polarization 
of air/water interface. With the symmetrical polariza- 
tion curves as shown in Fig^ such small signal can be 
determined with much better accuracy than previously 
reported results, where input polarization were only var- 
ied between 0° to 90°i2ii£ The polarization data with full 
360° scan range also provides another advantage. It is 
that the offset of the polarization angle in the SHG ex- 
periment, which can cause significant error on the very 
small SHG signal in sp polarization, can be accurately 
determined. 
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FIG. 1: 360° polarization SHG measurement of the neat 
air/water interface at two fixed detection polarizations (open 
circle for p detection; closed circle for s detection). Each point 
was averaged for three times. The lines are fitting curves. 



The data in Fig^ should follow the following 
functions^ 



I p = [A cos 2 9 + f?sin 2 9] 2 

I s = [Csm{29)} 2 (15) 

Fitting results for thep detection curve gave A = 14.42 ± 
0.07, B = 2.42 ± 0.18, the offset of 9 is 9 = 6.6° ± 0.4°; 
and the s detection curve gave C = 7.17 ± 0.08 and 
6» = 7.2° ± 0.5°. Such fitting results gave the follow- 
ing values, I pp — 208 ± 2, I sp — 5.9 ± 0.9 and /450s = 
51.4 ± 1.1. With the s and p detection ratio of 1.31, and 
after corrections with the Fresnel factors and the local 
field factors (which is essentially no influence), we obtain 
Xzxx/Xxzx = 0.30 ±0.03 and Xzzz/Xxzx = 2.56 ± 0.06. 
These ratios are consistent with the results reported in 
the literatures (Table HJ),and the values obtained here 
should be with better accuracy, as we compared the qual- 
ity of SHG data with those in the literatures, especially 
for the very small value of I sp £i£iiLH Nevertheless, be- 
cause I sp is very small, its value might have been overes- 
timated in the past literatures as discussed above. 

It is to be noted that in calculating Xzxx/Xxzx and 
Xzzz/Xxzx ratios, the assumption that Xzxx, Xxzx and 
Xzzz all have the same signs was used as in the previ- 
ous literatures Goh et al. studied these relative 
phases^i and our results are consistent with their re- 
port on this issue. 

B. Water Molecular Hyperpolarization Tensors 

Here we try to estimate the water molecular orienta- 
tion at the neat air/water interface from the SHG data. 

In Ea llOl and Ea llll there are three unknown param- 
eters, i.e. r, s and D. In order to get the value of D, 
value of either r, s or r/s needs to be known a priori. 
Otherwise, D value can only be estimated for possible 
r/s values. 

As discussed in Section II. C, (3 CCC term can be very 
small or negligible among the 7 nonzero terms of the Cz v 
water molecule. Now we plug in Xzxx/Xxzx = 0.30 ±0.03 
and Xzzz / Xxzx = 2.56 ± 0.06 in Tabled into EaHHand 
EcilTTfl and we obtain R = 0.69 ± 0.02 and D = 1.72 ± 
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0.02 by solving these two equations. If a 5-distribution 
function of 9 is assumed, we have 9 = arccos[(l/D) 1//2 ] = 
40.3° ± 0.4°. The calculation results with the previously 
reported Xzxx/Xxzx and Xzzz/Xxzx values are also listed 
in Tabled The R and D values obtained from previously 
reported data are all quantitatively consistent with our 
results. 

The result R < 1 confirms our expectation that (3 aac — 
fiaca or Pbbc = ftbcb are dominant terms, and the Klein- 
man symmetry is broken. The fact, that Xzxx/Xxzx is 
significantly smaller than R, indicated that the degree 
of broken of the Kleinman symmetry is usually not the 
same for the microscopic polarizability tensors and the 
macroscopic susceptibility tensors. From Franken and 
Wardj2& R value for a certain fundamental frequency (u) 
smaller than the resonant frequency can be estimated 
with the known low-lying electronic transition frequency 
(w e9 ). Franken and Ward defined^ 5 . 

I (2) _ (2), 

\A-ijk *-jik\ c Xijk 

-4^ ~ 1 (16) 

Clearly, R = 1/(1 + e). In our calculation, with laser 
wavelength as 800nm, and water transition wavelength as 
124nm, we obtained R = 0.86; with laser wavelength as 
800nm, and water transition wavelength as 190nm, R = 

0. 70. These R values well agree with our experimental 
value. 

Theoretical calculations have generated scattered val- 
ues for ratios of hyperpolarizability tensors for the wa- 
ter molecule, when the possibility for broken of Klein- 
man symmetry was considered £££±£££2^ J n many of 
other calculations, Kleinman symmetry was used as a 
constraining condition. [3 CCC value in all these calcula- 
tions are usually comparable or even larger than other 
tensor elements. However, results vary significantly with 
calculation models usedd So far, experimental values 
was scarce for comparison with these calculations No 
experimental value for ratio of hyperpolarizability ten- 
sors has been available. The experimental value of R 
we obtained here can be used for such comparison with 
calculated results. 

In our calculation, if [3 CCC can not be neglected, r, s and 
D can not be uniquely solved with Eq^Jand Eq^] If we 
let r = 0.5, which means 2f3 ccc = (3 caa + (3 c bb, then we get 
s = 0.95±0.04, i.e. R = 0.53±0.02, and D = 1.81±0.02, 

1. e. 9 = 42.0° ± 0.3°, from Xzxx/Xxzx = 0.30 ± 0.03 and 
Xzzz/Xxzx = 2.56 ± 0.06. If we allow (3 CCC to be even 
larger, for example, we let r = 0.25, then we get s = 
0.58±0.03, i.e. R = 0.43±0.02, and D = 2.08±0.03, i.e. 
9 = 46.1°±0.4°. These values indicated that even though 
r, s and D values can not be uniquely determined with 
experimental values of Xzxx/Xxzx and Xzzz/Xxzx, the R 
and D values are not very sensitive to the nonzero (3 CCC 
values used. Especially, the calculated 9 value changed 
only less than 6 degrees. 



Here we discussed how to extract polarizability tensor 
ratio R, or r and s, from the SHG experimental data con- 
sidering of dipolar (local or interface) contribution only. 
The values we obtained are self-consistent, and can pro- 
vide a good microscopic description of the SHG polar- 
ization signal from the neat air/water interface. These 
analyses indicate that quadrupole (nonlocal or bulk) con- 
tributions to SHG signal from air/ water interface can 
be very small or negligible, even though we can not say 
that such contributions has to be negligible. Most im- 
portantly, these analyses demonstrate that the previous 
rationales and evidences^**!*^ can no longer be used to 
draw conclusion for significant contributions of the bulk 
terms to the SHG signal from the neat air/ water inter- 
face. 



C. Water Molecule Orientation at Neat Air /Water 
Interface 

Here we discuss the orientation of water molecule at 
the neat air/water interface in the light of the SHG anal- 
ysis above. 

As we have known, there are at least two kinds of wa- 
ter molecules at the neat air/water interface from the 
SFG-VS experimental studies. 26 ! 65 ! 66 ! 67 ! 68 Namely, one 
kind (straddle-type) is with one free O-H bond stick- 
ing out of the liquid interface and with its dipole vec- 
tor (from positive to negative) point almost parallel to 
the air/water interface; and the other kind (non-staddle- 
type) is with both O-H bond hydrogen-bonded and the 
dipole vector point away from the liquid phasei££i§£ For 
straddle-type interfacial water molecule, both the free 
O-H bond and the hydrogen-bonded O-H bond can only 
be treated as Coq V symmetry vibrationallyiSSLS While 
the non-staddle-type water molecule belongs to the C?. v 
symmetry 67,68 The vibrational spectra of the stretching 
modes of both types of interfacial water molecules have 
been identified and the molecular orientation has been 
studiedi 26 . 65 i 66 i 67 i 68 

However, it is that the electronic transitions are re- 
sponsible for the SHG responses. Electronically, it is ex- 
pected that the straddle-type water molecule should still 
possesses most signature of C2 V symmetry, because the 
hydrogen bonding to one of the O-H bond perturbs much 
less significantly to the electronic transition, which has 
much higher transition energy and much broader spec- 
tra, than to the vibrational transition. Therefore, it is 
easy to conclude that the straddle-type water molecule 
should give very small SHG signal, with the dipole vec- 
tor almost parallel to the interface. Du et al. have shown 
that the straddle-type water molecule consists about 20% 
of the interfacial water molecules at the neat air/water 
interfaced Thus, it is easy to see that the SHG sig- 
nal comes almost all from the non-straddle-type water 
molecule. So, the D values determined in Section IV. B 
above provide orientational information only for the non- 
straddle- type water molecule at the neat air/ water inter- 
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face. 

The orientational parameter D obtained above indi- 
cated that the dipole vector of the non-straddle-type in- 
terfacial water molecule oriented around an apparent ori- 
entational angle of 40° from the interface normal. How- 
ever, this angle was calculated assuming a S distribution 
function of the tilt angle 9. As Simpson and Rowlen 
have shown, such an orientational angle is close to the 
'Magic Angle' 9 = 39.2°p£°. Therefore, it may come from 
a relatively broad orientational distribution. Indeed, be- 
cause there are different kinds of non-straddle-type wa- 
ter molecule at the interface region, as clearly indicated 
by the broad vibrational spectra between 3100cto _1 and 
3500cm- 1 in SFG-VS studies, 26 i 65 i 66 t 67 i 68 a relatively 
broad orientational distribution for the non-straddle- type 
hydrogen-bonded water molecule is generally expected. 
However, the distribution width should not be extremely 
broad. Otherwise, SHG or SFG-VS signal from the 
air/water interface should have been too weak to be de- 
tectable. On the other hand, recent detailed SFG-VS 
analysis has demonstrated that the orientational distri- 
bution of straddle-type water molecule appears to be as 
narrower as less than 15° 

On the other hand, broad orientational distribution 
has been predicted for both the straddle-type and non- 
straddle-type interfacial water molecule from Molecu- 
lar Dynamics (MD) and Monte Carlo (MC) simulation 
studiesi 71 i 72 i 73 i 74 i 75 i 76 i 77 i 78 i 79 i 80 i 81 i 82 Sokhan and Tildes- 
ley discussed the apparent disagreement of the early 
SHG and SFG results for the neat air/ water interface 
structured Goh et al.'s early SHG studies indicated that 
water hydrogens are pointing towards the liquid side. 3 
While Du et al.'s SFG-VS study identified the staddle- 
type water molecule at the interface. From their simu- 
lation results, Sokhan and Tildesley believed that SHG 
and SFG-VS techniques arc sensitive to different parts of 
the interface, respectively. Simulation results often found 
that the plane of water molecule is aligned approximately 
parallel to the interface on the liquid side and perpendic- 
ular to the interface on the gas side, and the dominant 
contribution to the interface susceptibility is from the 
liquid side, 72 ! 73 ! 74 ! 75 ! 76 ! 77 ! 78 ! 79 ! 80 ! 81 ! 82 Unfortunately, such 
conclusions are drastically in disagreement with known 
SHG and SFG-VS experimental observations of the neat 
air /water interface, because when the water molecule has 
its molecular plane lying parallel to the interface, it is 
impossible to detect any SHG or SFG-VS signal in any 
polarization and experimental configuration. Detailed 
analysis of the SHG and SFG-VS experimental data has 
been significantly improved and systematically demon- 
strated in recent years, 43 ! 44 ! 55 ! 67 ! 68 ! 83 ! 84 ! 85 ! 86 ! 87 ! 88 ! 89 ! 90 ! 91 
Even though simulation results have provided a great 
deal of detailed understandings of the liquid interfaces, 
these disagreements between experimental and theoreti- 
cal results suggest that the physical picture provided by 
the simulations may need to be systematically reexam- 
ined as discussed below. 

As we have demonstrated above, SHG probes mainly 



the non-straddle- type of interfacial water molecule, which 
has both O-H bond hydrogen-bonded. On the other 
hand, SFG-VS probes the stretching vibrational spectra 
of both the straddle-type and the non-straddle-type in- 
terfacial molecules. 26 - 65 - 66,67 - 68 Recent detailed analysis 
of the SFG-VS spectra in different polarizations and ex- 
perimental configurations have helped to clarify the spec- 
tral assignments, symmetry determination of the interfa- 
cial water species, as well as understanding more on the 
orientation and motion of water molecules at the neat 
air/water interface. 6 ' -68 It has been known that most 
of these non-straddle-type water molecules have either 
4 or 3 hydrogen-bonds on average .Si&SLSiSSi Therefore, 
these non-straddle-type water molecules have to reside 
on the liquid side instead of the gas side as suggested by 
the simulation results. As probed by SHG and SFG-VS, 
their molecular plane can not be close to parallel to the 
interface. Of course, it is clear that the straddle-type 
water molecule is at the outermost edge of the liquid 
boundary, and its molecular plane is certainly close to 
perpendicular to the liquid interface. However, to call 
this layer as the gas side of the interface may not be ap- 
propriate, because the air/water interface is extremely 
sharp (3.2A) from X-ray reflectivity 9 ^*^ and theoreti- 
cal simulationsdd As the experimental analyses have 
consistently given relatively narrower orientational dis- 
tribution of the liquid interfacial molecules, simulation 
results have always accounted for a much broader orien- 
tational distribution, 82 ! 86 ! 87 ! 89 ! 90 ! 91 ! 95 ! 96 ! 97 Such issue has 
to be answered in the future studies. 

Here we make one final comment on the orientational 
energy of the air/water interface. Using the temperature 
effect of SHG signal they measured, Goh et al. concluded 
that it is about \kT, which suggested a quite disordered 
air/water interface at the room temperature,— Since now 
the temperature effect has to be gone, the value of the 
orientational energy from the same treatment of Goh et 
al. has to be much bigger than a few kT. This neverthe- 
less suggests a much better ordered air/ water interface, 
and it is consistent with the hydrogen-bonding nature of 
the molecules at the neat air /water interfaced The hy- 
drogen bond energy is typically about 5kcal/mol, and it 
implies a much well ordered interface structure. 

In short conclusion, non-resonant SHG provides ori- 
entational measurement of the non-straddle-type of in- 
terfacial water molecule. Its dipole vector orientation is 
distributed around 40° from the interface normal. This 
picture is fully consistent with and complementary to 
the SFG-VS experimental results. Molecular simulation 
results may need to be reexamined for detailed orienta- 
tional structure of the neat air/water interface. 



V. CONCLUSION 

As Shen had insightfully pointed out, the major ob- 
stacle in the advance of the (SHG) technique for surface 
studies is the lack of microscopic understanding of the 
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surface nonlinear susceptibility Xs~ However, the over 
simplifying uniaxial molecular model has been gener- 
ally used in the field. In recent years, efforts for sys- 
tematic treatment of the problem beyond the uniaxial 
molecular model have been attempted by Simpson and 
co-workers i Our recent experiences in quantitative 
treatment of the SFG-VS spectra, which always involves 
treatment of non-uniaxial molecular groups, prompted us 
to work similarly on the treatment of SHG of molecular 
interface.: 44 - 55 - 67 - 68,83 - 84 - 85 - 86 - 87 - 89 - 90 - 91 

In this work, we re-examined the problem of the 
SHG reflected from neat air/water interface with detailed 
treatment on the connection between the microscopic po- 
larizability of water molecule and macroscopic suscepti- 
bility of the interface. Because these treatment on mi- 
croscopic polarizability is based on the molecular sym- 
metry properties of single water molecule, so it can be 
well suitable for treating water molecule, as well as other 
molecules with the same symmetry, at different molecu- 
lar interfaces. After examination on issues of the origin 
and validity of the Kleinman symmetry, we have con- 
cluded that the broken of macroscopic Kleinman sym- 
metry of the neat air/ water interface is not an indication 
for quadrupole (bulk) contribution to the surface SHG 
signal, as previous studies suggested and generally be- 
lieved in the field. We further demonstrated that using 
dipole contribution only can fully address the broken of 
macroscopic Kleinman symmetry of the air/ water inter- 
face. Using properly combined molecular polarizability 



tensor terms, we further simplified the treatment and 
are able to obtain the value of orientational parameter D 
from SHG data. Such procedures can be generally trans- 
ferred to the treatment of interfacial molecules belong to 
different molecular symmetry. 

In conclusion, the SHG signal from the neat air/water 
interface can be treated fully with dipole (local or inter- 
face) contribution. The broken of macroscopic Kleinman 
symmetry is not a sufficient condition for quadrupole 
(nonlocal or bulk) contributions in surface SHG. Even 
though in general the quadrupole (bulk) contributions 
has to be considered in SHG studies, this work sug- 
gests that they can be very small or even negligible for 
the neat air/water interface, which has long been con- 
sidered the primary case favoring dominant quadrupole 
(bulk) contributions. The successful treatment of the 
neat air/water interface with dipole contribution terms 
indicates that SHG probably is indeed a probe of the 
surfaces and interfaces of isotropic fluids, as strongly ar- 
gued by Andrews et a ;^iLigii2iS& The treatment and 
approaches we presented here can be applied and tested 
with other fluid interfaces. 
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